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Analytical expressions of the one- and two- body terms in the cluster expansion of the one-body 
density matrix and momentum distribution of the s-p and s-d shell nuclei with N = Z are derived. 
They depend on the harmonic oscillator parameter b and the parameter /3 which originates from 
the Jastrow correlation function. These parameters have been determined by least squares fit to 
the experimental charge form factors. The inclusion of short-range correlations increases the high 
momentum component of the momentum distribution, n(k) for all nuclei we have considered while 
there is an A dependence of n(k) both at small values of k and the high momentum component. 
The A dependence of the high momentum component of n(k) becomes quite small when the nuclei 
24 Mg, 28 Si and 32 S are treated as ld-2s shell nuclei having the occupation probability of the 2s-state 
as an extra free parameter in the fit to the form factors. 
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The momentum distribution (MD) is of interest in many research subjects of modern physics, including those 
referring to helium, electronic, nuclear, and quark systems In the last two decades, there has been significant 

effort for the determination of the MD in nuclear matter and finite nucleon systems @-^] . MD is related to the cross 
sections of various kinds of nuclear reactions. Specifically, the interaction of particles with nuclei at high energies, 
' such as (p,2p), (e,e'p), and (e,e') reactions, the nuclear photo-effect, meson absorption by nuclei, the inclusive proton 
&\ , production in proton-nucleus collisions, and even phenomena at low energies such as giant multipole resonances, 

■ give significant information about the nucleon MD. The experimental evidence obtained from inclusive and exclusive 
electron scattering on nuclei established the existence of a high- momentum component for momenta k > 2 fm _1 
pjj p7{ . It has been shown that, in principle, mean field theories can not describe correctly MD and density distribution 
simultaneously JlJ] and the main features of MD depend little on the effective mean field considered |D| . The reason 
is that MD is sensitive to short-range and tensor nucleon-nucleon correlations which are not included in the mean field 
theories. Thus, theoretical approaches, which take into account short range correlations (SRC) due to the character 

'"pi \ of the nucleon-nucleon forces at small distances, are necessary to be developed. 

Zabolitzky and Ey jlj, employing the coupled-cluster (or exp(S')) method for the microscopic evaluation of nuclear 
MD for the ground states of 4 He and 16 and using various realistic NN-potentials, showed that the contribution of 
correlations dominates for momenta beyond 2 fra — 1 . A realistic interaction and a many-body approach have been 
used by Benhar et al |[3) for the evaluation of MD of 12 C, 16 O and 40 Ca. Their results have yielded a much larger 
. content of the high momentum component with respect to the results obtained within the Hartree-Fock approach or 
J^j ' within methods which take into account the effect of correlations phenomenologically. 

Bohigas and Stringari || and Dal Ri et al Q evaluated the effect of SRC's on the one- and two- body densities by 

■ developing a low order approximation (LOA) in the framework of Jastrow formalism. They showed that one-body 
quantities provide an adequate test for the presence of SRC's in nuclei, which indicates that the independent-particle 
wave functions cannot reproduce simultaneously the form factor and the MD of a correlated system and also the effect 
of SRC's strongly modify the MD by introducing an important contribution in the region k > 2 fm _1 . Stoitsov et al 
fl7[ generalised the model of Jastrow correlations within the LOA of Ref. pq| , to heavier nuclei as 16 0, 36 Ar, 40 Ca. 
Their analytical expressions for the MD show the high momentum tail. They found that there is an A dependence 
of MD for small values of k, while for large values of k the slope of logn(fc) versus k is roughly the same for the 
above three nuclei as well as for 4 He. The same behaviour of the MD of protons and neutrons for the A = 3, 4 nuclei 
has been found earlier by Schiavilla et al jl^] performing variational calculations with realistic interactions. MD for 
the nuclei 4 He, 16 and 40 Ca was also calculated by Traini and Orlandini jll} within a phenomenological model in 
which dynamical short-range and tensor correlations effects were included. They showed that SRC increase the high 
momentum component considerably while the tensor correlations do not affect the MD appreciably ]ll],p9[- I n heavy 
nuclei, the local density approximation was used |L6| for the study of the effect of SRC's in MD and the predictions 
were in agreement with the results of microscopic calculations in nuclear matter and in light nuclei. 

The influence of SRC's on the MD of nucleons in nuclei has also been evaluated by Muther et al pl| within the 
Green-function approach assuming a realistic meson-exchange potential for the nucleon-nucleon interaction. Their 
analysis on 16 O demonstrates that a non-negligible contribution to the MD should be found in partial waves which 
are unoccupied in the simple shell model. Another approach is to consider the average occupancy of the relevant 
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shell-model orbitals ||. It has been found that the depletion of such orbitals can be of the order 15% or more for 
single-particles (SP) states below the Fermi energy fl30f . 

In the various approaches, the MD of the closed shell nuclei 4 He, 16 and 40 Ca as well as of 208 Pb and nuclear 
matter is usually studied. There is no systematic study of the one body density matrix (OBDM) and MD which 
include both the case of closed and open shell nuclei. This would be helpful in the calculations of the overlap integrals 
and reactions in that region of nuclei if one wants to go beyond the mean field theories J3^] . For that reason, in the 
present work, we attempt to find some general expressions for the OBDM p(r,r') and MD n(k) which could be used 
both for closed and open shell nuclei. This work is a continuation of our previous study pl| on the form factors and 
densities of the s-p and s-d shell nuclei. The expression of p(r, r') was found, first, using the factor cluster expansion 
of Clark and co-workers [[35-37 and Jastrow correlation function which introduces SRC for closed shell nuclei and 



then was extrapolated to the case of N = Z open shell nuclei. n(k) was found by Fourier transform of p(r,r'). 
These expressions are functionals of the harmonic oscillator (HO) orbitals and depend on the HO parameter b and 
the correlation parameter (3. The values of the parameters b and /3, which we have used for the closed shell nuclei 
4 He, 16 and 40 Ca, are the ones which have been determined in Ref. |3l|] by fit of the theoretical F c h{q) 7 derived with 
the same cluster expansion, to the experimental one. For the open shell nuclei 12 C, 24 Mg, 28 Si and 32 S we provide 
new values for these parameters, which have been found to give a better fit to the experimental form factors than in 
our previous analysis jpj. It is found that the high- momentum tail of the MD of all the nuclei we have considered 
appears for k > 2 fm~ and also there is an A dependence of the values of n(k) for 2 fin - < k < 5 fm _ . This A 
dependence of MD was first investigated considering 24 Mg, 28 Si and 32 S as Id shell nuclei. Next we treated the above 
nuclei as ld-2s shell nuclei having the occupation probability of the 2s state as an extra free parameter in the fit of 
the form factors. The A dependence is quite small in the second case. 

The paper is organised as follows. In Sec. II the general expressions of the correlated OBDM and MD are derived 
using a Jastrow correlation function. In Sec. Ill the analytical expressions of the above quantities for the s-p and s-d 
shell nuclei, in the case of the HO orbitals, are given. Numerical results are reported and discussed in Sec. IV, while 
the summary of the present work is given in Sec. V. 

II. CORRELATED ONE-BODY DENSITY MATRIX AND MOMENTUM DISTRIBUTION 

A nucleus with A nucleons is described by the wave function Vf^ri, r 2 , r^) which depends on 3 A coordinates as 
well as on spins and isospins. The evaluation of the single particle characteristics of the system needs the one-body 
density matrix |33[|3~i| ] 

p(r,r') = yV(r,r 2 ,...,r A ) *(r', r 2 , • • • , r A ) dr 2 • • ■ dr A , (1) 

where the integration is carried out over the radius vectors r 2 , • • • , ta and summation over spin and isospin variables 
is implied. p(r,r') can also be represented by the form 

p(r,r') = -^|p = N{9\0 Tr ,\V) = N(0„,) , (2) 

where VP' = , i , (r' 1 ,r 2 , ■■■,r' A ) and N is the normalization factor. The one-body "density operator" O rr /, has the form 

A A 

O rr , = ]T 5( ri - r)5(r> - r') J[ S(r 3 - r$) . (3) 

i=l jjti 

In the case where the nuclear wave function 'J can be expressed as a Slater determinant depending on the SP wave 
functions 4>i(r) we have 



PSD 



(r,r') = $>*(r)<Mr') (4) 



The diagonal elements of the OBDM give the density distribution 

p(r,r) = p(r) , (5) 
while the MD is given by the Fourier transform of p(r,r'), 
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i(k) = J exp[ik(r- r')] p(r,r') dr dr' . 



(6) 



In the case of a Slater determinant, MD takes the form 

A 



n SD (k)=^0*(k)^(k) , (7) 



i=i 

where 



4(k) = — ^ J &(r)exp[ikr]dr. (8) 
The second moment of the MD is related to the expectation value of the kinetic energy, (T), by the expression 

(T) = [ n(k)fc 2 dk. (9) 
2m J 

A. One-body density matrix 

If we denote the model operator, which introduces SRC, by J 7 , an eigenstate $ of the model system corresponds to 
an eigenstate 

* = F§ (10) 

of the true system. 

Several restrictions can be made on the model operator as for example, that it depends on (the spins, isospins 
and) relative coordinates and momenta of the particles in the system, that be a scalar with respect to rotations etc. 
PH . Further, it is required that T be translationally invariant and symmetrical in its arguments \- A and 

possesses the cluster property. That is if any subset i\ ■ ■ ■ i p of the particles is removed far from the rest i p +i ■ ■ ■ ia, F 



decomposes into a product of two factors, J- (I ■ ■ ■ A) = T{i\ ■ ■ -i p ) T(i p +i ■ ■ ■ ia) |37|. In the present work T is taken 
to be of the Jastrow type (3j| , 

A 

^ = U.f(ra), (11) 

i<j 

where f(rij) is the state-independent correlation function of the form 

/(r y ) = 1 - exp[-/3(r, - r,-) 2 ] . (12) 

The correlation function /(ry) goes to 1 for large values of =| — rj | and it goes to for — ► 0. It is obvious 
that the effect of SRC, introduced by the function /(rjj), becomes large when the SRC parameter (3 becomes small 
and vice versa. 

In order to evaluate the correlated one-body density matrix p eor (r,r'), we consider, first, the generalized integral 

J(a) = (*|exp[a/(0)O Pr /]|* / ) , (13) 
corresponding to the one-body "density operator" O rr / (given by (^)), from which we have 

~d In 1(a)" 



(Orr,) = 



da 



(14) 



For the cluster analysis of equation (|14|), we consider the sub-product integrals |35|-|37|, for the sub-systems of the 
A-nucleons system 
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h(a) = (i | ^ t (r 1 )exp[a/ l (0)o rr ,(l)]^(r' 1 ) | i') , 
Uiia) = (ij | ^ t (r 12 )exp[a^(0)o rl .(2)]^(r' 12 ) | i'j'), , 

Iijk{oc) = (tjfc | ^ t (r 12 )^t( ri3 )^t (r23)exp [ a ^ fe(0)Orr , ( 3 ) ]^ ( / i2) ^ (r ' i3) ^ ( ^ 3 ) | ,y fc ') £ 



Ju-.a = 1(a) , (15) 
where the operators o rr /(l), o rr /(2), . . . have the form 
o rr /(l)=5(n-r)<5(ri-r') , 

o rr ,(2) = <J(r x - r)*(rj - r')S(r 2 - r' 2 ) + S(r 2 - r)S(v' 2 - r')5(n - r' x ) , (16) 

and so on. The factor cluster decomposition of the above integrals, following the factor cluster expansion of Ristig,Ter 
Low, and Clark p5|-E$7) , gives 



<O pp ») = <0„/)i + (O rr ,) 2 + • ■ ■ + (O rr ,) A , (17) 

where 

A r . 



<o rr <)i = E 



din Ii(a) 
da 



(18) 



A „ 

(O rr ,) 2 = — pn/ i3 -(a) - In/i(a) - \nI 3 (a)] a=Q , (19) 

and so on. T(r\) is chosen to be the identity operator. 

Three- and many-body terms will be neglected in the present analysis. Thus, in the two-body approximation, 
p cor {r,r'), defined by Eq. (||), is written 

p cor (r, r') » iV[(O rr ,)i + (O rr ,) 22 - (0„/) 3 i] , (20) 

where 

A 

(O rr /)i = ^(i | o rr /(l) | i') , (21) 



O rr ,} 22 = | ^( ri2 ) 0rr ,( 2 )^( r ; 2 ) | f/) a| (22) 



(O rr ,) 2 i =2<« I °rr'(2) I (23) 

i<j 

If the two-body operator T(r' l2 ) is taken to be the correlation function given by Eq. (|lj), then 

T\r 12 )o rr ,(2)T{r' l2 ) = o„,(2) [1 - gl (r,r 2 ) - g 2 (r',r 2 ) + g 3 (r,r',r 2 )] , (24) 

where 

gi(r,r 2 ) = exp[-/3(r 2 + rf )] exp[2$rr 2 ], g 2 (r',r 2 ) =gi(r',r 2 ), 
g 3 (r, r', r 2 ) = exp[-/3(r 2 + r' 2 )] exp[-2/?r§] exp[2/3(r + r')r 2 ] (25) 

and the term (O rr /) 22 is written 
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(O rr ') 2 2 = (O rr ') 2 i - 22 (r,r',gi) - 22 (r,r',g 2 ) + 22 (r,r',g 3 



(26) 



where 



C»22(r,r',g,) = J^iij I o rr ,(2)g,(r,r',r 2 ) | i'f) 



= 1,2,3 



i<3 



Thus, p cor (r,r') takes the form 

p cor (r, r') a iV[(O rr -)i - 22 (r, r', gl ) - 22 (r, r', g 2 ) + 22 (r, r', g 3 )] 
This is also expressed in the following form 

p cor (r,r') AT pgr, (r, r') + / [-gi(r, r 2 ) - g 2 (r', r 2 ) + g 3 (r, r', r 2 )] 
x [p S D(r,r')p SD (r 2 ,r 2 ) - p SD (r, r 2 )p SD (r 2 , r')] dr 2 



(27) 



(28) 



(29) 



where psi>( r , r ') is the uncorrelated OBDM associated with the Slater determinant. 

It should be noted that a similar expression for p cor (r,r'), given by Eq. (29), was derived by Gaudin et al. [^8| in 
the framework of LOA. Their expansion contains one- and two-body terms and a part of the three-body term which 
was chosen so that the normalization of the wave function was preserved. Expression (|2^) of the present work has 
only one- and two-body terms and the normalization of the wave function is preserved by the normalization factor N. 

In the above expression of p cor (r,r'), the one-body contribution to the OBDM is well known and is given by the 
equation 



(O rr ')i = p SD (r,r') = ~y^r] n i(2l + l)^* ; (r)0 n ; (r')Pi (cos u rr >) 



(30) 



where r) n i are the occupation probabilities of the states nl (0 or 1 in the case of closed shell nuclei) and 4> n i (r) is the 
radial part of the SP wave function and oj rr i the angle between the vectors r and r'. 

The term 22 (r, r', gi), performing the spin-isospin summation and the angular integration, takes the general form 

22 (r,r',g £ ) = 4 VmkVn.lA^i + 1)^ + 1) 



where 



^S^f (r,r',g.) - £ (kO^O I fc0) 2 <^; fc (r,r',g,) 



k=0 



, £=1,2,3, 



(31) 



KllTntlf ( r > r '>Si) = TZ&rhtr) tn 3h (r') exp[-/?r 2 ] P i3 (cosuw) 



■in 



x / K 2 i 2 ( r 2) tntufa) exp[-/3r|] i k (2(3rr 2 ) r\ dr 2 , 



and the matrix element A™ 3 i \ 3 i ™* l l *' k (r , r', g 2 ) can be found from ( |32"| ) replacing r < — > r' and nylx 
matrix element corresponding to the factor g 3 is 

(r,r',g 3 ) = ^K lh (r) exp[-/3(r 2 + r' 2 )] ^(aw) x 



K-ufo) </Vz 4 (r 2 ) exp[-2/3r|] i k (20\r + r'|r 2 ) r\ dr 2 , 



(32) 

77,3/3 while the 



(33) 



In Eqs. ( |32| ) and ( |33| ) the modified spherical Bessel function, i k (z), comes from the expansion of the exponential 
function exp[2/3xix 2 ] of the factors g£ in spherical harmonics, that is 



exp[2/3 Xl x 2 ] = 4vr^ i k (20xix 2 )Y£ mk (Q 1 )Y kmk (Sl 2 ), 



km 1, 
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while the factor Clf t (u> rr i) which depends on the directions of r and r' is, 



^hhi^rr') — ^2 



(h - mi/ 3 m 3 | kO) ( (I i - mi)! (l 3 - m 3 )l 



W I feO) l Px+rm)ia, + m,)l J P <> (C ° S ^ } P '« (C0S ^' } ' 



(34) 



where w r and w r < are the angles between the vectors r, r + r' and r', r + r', respectively. The final expression of 
fif iZ3 (w rr <) depends on cu rr >, the angle between the vectors r and r'. 

The expression of the term 022(1", v',ge) depends on the SP wave functions and so it is suitable to be used for 
analytical calculations with the HO orbitals and in principle for numerical calculations with more realistic SP orbitals. 
Expressions ( p0| ) and ( J3l| ) were derived for the closed shell nuclei with N = Z, where rj n i is or 1. For the open shell 
nuclei (with N = Z) we use the same expressions, where now < r\ n i < 1. In this way the mass dependence of the 
correlation parameter (3 and the OBDM or MD can be studied. 

Finally, using the known values of the Clebsch-Gordan coefficients, Eq. (|3l|), for the case of s-p and s-d shell nuclei, 
takes the form 



22 (r,r',g f ) = 

o 10000,0/ / \ 2 , Too ,(0101.0/ / \ r ,(0101.2/ / \1 2 , o /i 1010.0/ / 

3^oooo ( r , r ,gt)m.+ 33A 0lnl (r,r ,g £ ) -6A 0101 (r,r,g/) rj lp + 3A 1010 ' (r,r,g^7? : 



72s 



95A 



0202,0/ 
0202 



50 



90 



(r,r',g£) - — A 0202 > (r,r',g £ ) - — A 0202 ' (r,r',g 



1040001,0/ / s ,-,040100,0/ / x _ 040100,1/ / x _ o^OOOl.l/ / 
i/yi oooi l r ' r ! Se) + lz/i oioo V r ; r ;S^; d A)ooi \ ' , g^ j 0^0100 \ r > r i&e 



20A 



0002,0 
0002 



0200,0 
0200 



(r,r',g,)-5A°^ 2 (r,r',g 



5A 



0002,2 
0200 



(r,r',g £ 



VuVid 



0010,0, 
0010 



/ / \ , a j 1000,0/ I \ ,.1000.0/ / \ ,(0010,0/ / \ 

(r, r , gt) + 4A 1000 ' (r, r ,g e )- A 0010 ' (r, r ,g e )- A 1000 ' (r, r , g e ) 



an a 0102,0/ / \ , C r\ A 0201,0/ / \ r ^0201,1/ / \ r A 0W2A / 1 \ 
60 A)102 0, r , Et) + 60 A)201 ( r > r > St) ~ 6 A)102 ( r > r > S<0 ~ 6 A)201 ( r > r >gf) 



9A 



0201.3 
0102 1 



r,r',g,)-9A™(r,r',g. 



VlpVld 



-,0^0110,0/ / \ . -, n ,(1001,0/ / \ o A 1001,1 / / \ o A 0110,1 / / \ 

12^0110 (r,r,g/) + 12^00!' (r, r , g t ) - 3A ono ' (r, r , gt) - 3A W0 { (r,r ,ge) 



r\\ v r\2s 



20,4 



0210,0/ 
0210 



/ / \ 1 on /( 1002,0/ I \ r ,(1002,2/ / \ r ,(0210,2/ / \ 

(r, r , g e ) + 20A 1002 (r, r , g e ) - 5^ 0210 ' (r, r , gt) - 5A 1002 ' (r, r , g e ) 



^210,2/ 



(35) 



It should be noted that Eqs. (^Tj) and (pq) are also valid for the cluster expansion of the density distribution and 
the form factor as it has been found in ref. |H| and also in the cluster expansion of the MD. The only difference is 
the expressions of the matrix elements A. 



B. Momentum distribution 



The MD for the above mentioned nuclei can be found either by following the same cluster expansion or by taking 
the Fourier transform of p(r,r') given by (|2q). In both cases the correlated momentum distribution takes the form 



rw(k) (6 k )i - 20 2 2(k,gi)+0 2 2(k,g 3 ; 



(36) 



where 



(6 k ) 1 = n SD (k) = ^0j(k)^(k). 



(37) 



i=l 



G 



The term 22 (k, g^), as in the case of OBDM, is given again by the right-hand side of Eqs. ( |3l| ) and ( |35| ) replacing 
\n3i3n4i4,' 

nilm2h 



the matrix elements j4™ 3 j 3 ™ 4 ; 4 ' fc (r, r', g^), defined by Eqs. ([32]) and (|33|), by the Fourier transform of them, that is by 



the matrix elements 

AlfCif^gi) = (2^)3 / ^^^'(r.r'.g^exptikCr-rOldrdr', t= 1,3 . (38) 

As in the case of the OBDM, expression ( |36| ) is suitable for the study of the MD for the s-p and s-d shell nuclei 
and also for the study of the mass dependence of the kinetic energy of these nuclei. The mean value of the kinetic 
energy has the form 

(T) = N[(T) 1 - 2T 22 ( gl ) + T 22 (g 3 )] , (39) 

where 

(T)i = ^ J fc 2 n SD (k)dk , T 22 (g,) = A- J fc 2 22 (k, g,)dk , I = 1, 3 . (40) 



III. ANALYTICAL EXPRESSIONS 

In the case of the HO wave functions, with radial part in coordinate and momentum space, 

^(r) = N nl b-^r[L 1 ^ (rf) e^ 2 , r b = r/b , 

Mk) = i l (-l) n+l N nl b 3 / 2 k l b L 1 ^ (k 2 b ) e"^ 2 , h = kb , 

where 



(41) 



Nr, 



2n\ 



T(n + l 



1/2 



analytical expressions of the one-body terms, (O rr ')i an d (Ok)i as well as of the matrix elements A"^ 3 ™ 4 ^ 4 ' (r, r', g^) 
and A™** 1 ^*' (k, gt), which have been defined in Sec. II, can be found. From these expressions, the analytical 
expressions of the terms 22 (r, r',ge) and 22 (k, g^), defined by Eq. (|35|), can also be found. 
The expressions of the one-body terms, (O rr /)i and (Ok)i, have the forms 



(O rr ')i = p S D(r,r') 



tt 3 / 2 6 3 



2rj ls + 3r] 2s - 2ij 2s {r 2 + r' b ) + 4:T] lp r b r b cos u. 



+ g[»?2a + ?7id(3cos 2 uj rr , - l)r 2 r b 2 } 



exp[-(r 2 +^)/2] 



(42) 



where the coefficients C 2 fc are 



2b 3 



6k) x = wsc(k) = -^exp[-fc 2 ] ^C 2fc fc, 



2A- 
b i 



fe=0 



(43) 



C = 2r] U + 3r? 2s , C 2 = 4(7?i p - J7 2s ) , C 4 = -(2^ + r) 2s ) ■ 
The analytical expressions of the matrix element A™ 3 |^ 4 ' 4 ' fc (r, r', g^) [I = 1,3) have the form 



l + 3y 
2(1 + 2/)' 



< J}, (cos aw) ^ ^ ^ 2U)4 (y) LjJ +l4 _ fc 



U> 2 = W4—O 



-V 2 2 



k)+w 2 +w 4 I J _|_ „.' b I ' 



(44) 



(45) 
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and 



<Kf (r,r',g 3 ) = B b-*y k \v b + r' b \ k r<V b ' 3 ^'(rg) C 8 ^)"!^^) 



x exp 



exp 



l + 2y 



n 2 Hi 



W2— u;4— 



3((2+U— fe)+«j 2 +«J4 V 1 + 2?/ 



-(n + r' b ) 



(46) 



where ?/ = (3b 2 and 



Bn = 



(47) 



B W2Wi {z) = [l(h + k-k)+w 2 +w 4 }\ J] L^fa+^ + s) (! + ,)-*»*-.- 1^) 



(48) 



while the one corresponding to the factor g2 can be found from ( |45| ) replacing rf, < — > and niZi < — > ^3/3. 

The substitution of A™^™^* ,,c (r, r', g^) to the expression of 022(1", v'i9t) which is given by Eq. ( |35| ) leads to the 
analytical expression of the two-body term of the OBDM, which is of the form 



022(rb,r&) = /i(r & ,r£,cosav r ')exp 
+/iM,,r 6 ,cos w„.,)exp 
+h{n,r' b , cos w rr /)exp 



l + 3y 
"2(l + y) 5 

1 + 3?/ ,2 _ 1 2 
2(l + y) 7 ' h 2 r " 

1 + 22/ 2 , ,2, 



exp 



(r b + r^) 2 



l + 2y 



(49) 



where ft (r b , r' b , cos w rr ' ) , (•£ = 1,3) are polynomials of r b , r' b and coso-Vr' which depend also on y = /3b 2 and the 
occupation probabilities of the various states. 

The corresponding analytical expressions of the matrix elements -A" 3 ^ 3 ™ 4 ^ 4 ' (k, g,«), (£ = 1,3) which contribute to 
the two-body term of the MD were found substituting <fr n i (r) with that of the HO wave function into Eq. ( p8[) . The 
expression of -A™?;^ 4 ;^ (k, gi), which can be found easily, has the form 



nil 1 n 2 l2 



(k, gl ) = B Q b 3 (-l) n *k 2 b h L l £ h (k 2 b ) exp 



l+2y j 
1 + 3?/ 



(.V) 



uii— W2—O W4—O 



t=0 



-(h-h + k) + wi+t) 



( ( "1 + ^1 + \ \ 2 |(il-i3)+Wl n +y)3(il+i3-i2-U)+«>l-'U>2-U>4 



:(1 + 3y )-i(3+ii+?3+fc)--i-* L ^ 



(h-h+k)+ Wl +t \^2(1 + 3y) 



1+y -kl 



(50) 



where 



B^tG/) = (V2y) fe+2t 
1 



-(l + / 2 + ^4 + fc) + W2 + W4 



k — li — Id 



IL>2 — W4 



{k + t- 



(-l)«* +t fn t + h + k 



(51) 



The expression of A 



nil±n2h 



(k, g3) is more complicated. It has the the general form 
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A 



^S' fc (k,g 3 ) = ^5o6 3 exp 



k 2 
K b 



«2 "4 2^ 

E 

t=0 



E 



W2— 1L>4— 



l + 2y 

x(l + 2y)~^ 3+ ' 2+ ' 4+fc -'~ tU2 ~"' 4- */" 3 ; ' 3, ' c (fc( ) ) 



(y) 



(52) 



The general expression of the quantity ^fjf' is quite complicated. For that reason we calculated it for various 
cases which are needed for the s-p and s-d shell nuclei. The various cases and the corresponding expressions of 
I2fii' k (kb) are given bellow. 



1. Case 1: nih = 713/3 and k = 

[il/2] ni p+a2(p+a-T)l 1 -2pn 1 -o-n 1 -a-p p , >.| Q+T+m +m 

c:'»=eee ESSE E — 



nilplalwilios! ~ a+Wl+T+v \l + 2y) 



p=0 ct=0t=0 l/=0 a=0 p=0 u>i=0 io 3 =0 

^1 - 2p\ ( p + a\ ( 2(p + a - r)\ / ni + a — p \ ( p, + lx — \ 



a I V T / V v J \ Ul ~ a ~ P ~ Wl J \ L 1 ~ w 3 

(a + w 1 +T + v)\ {2h - 2p)!(ni + h + + w 3 + t + 2cr - a - r - y + |)! 

X 2r + l G 1 -p)!(Z 1 -2p)!( ( T + Z 1 + i)! 



(53) 



r01,2 
'01 



£. Case 2: nih = 713/3 — 01 and k = 2 
(fc b ) = 2 1+ * (J + tj !(1 + 22/) 2+t (l + 4y)"i 



(54) 



3. Case 3: mil = n 3 / 3 = 02 and k = 2,4 



T 02,k, 
l 02 



(k b ) = 22 +t (l + 2y) L+t (l+4y)-2- t 

(i±i + 1 ) !( i + + %r t - (5 + «) IL? (ji- 



(55) 



4- Case 4-' nih 7^ ^3/3 and h = or/and I3 = 
"i "3 [ ; i/2] ['3/2] h-2p! l 3 -2p 3 pi+wi P3+W3 pi+p3+wi+w 3 -ai-a 3 



Qt k (h) - EEEEEEEE E 

wi- W3- pi— P3— Ti= r 3 — cri— 0-3— 



y=0 



n (-i) pi+u " 



+T3+(T 3 



2(n + r 3 + cti + 03 + 1) 



i=l,3 



- Pi)!^ - 2^)! V 



1 + fc + h + l 3 - ai - C7 3 - ri - 



+ t + wi + w 3 - v 



<Tl + CT 3 + Tl + T 3 



. Pi + p3 + + w 3 - cri - er 3 \ 

! 1/ J i(cr 1 +<J3+ri+T3+2iy) ^ _|_ 2y J 



(56) 
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5. Case 5: mh = 01, 713/3 = 02 or nih = 02, 713/3 = 01, and k = 1, 3 



C k (h) = C k (h) = 2H^i + | + ^ !( i + 2y )i(3+^(i + % ) 



-3-#-t 



„ k 

2 + -+t 



l + 4y 



k 2 + k-6 3J(l + 2y) 2 1 \l + 2y 



k 2 



(57) 



The substitution of -A™^^ 4 ; 4 '*^, gi) to the expression of (322(k, g^) which is given by Eq. ( |35| ) leads to the analytical 
expression of the two-body term of the MD, which is of the form 



22 (fc) = /i(fc b 2 )exp 



1 + 2y a 



/ 3 (fc 2 )exp 



l + 2y 



(58) 



where fi(k 2 ), (£ — 1,3) are polynomials of fc 2 which depend also on y = (3b 2 and the occupation probabilities of the 
various states. Similar expressions have been found for the moan value of the kinetic energy 

It should be noted that, although the above expressions of the matrix elements A and A seem to be quite complicated, 
they can easily be used for analytical calculations with programs such as Macsyma or Mathematica. As the above 
expressions have been found for the (N = Z) s-p and s-d shell nuclei they can be used for the systematic study of the 
OBDM and MD in this region of nuclei. 



IV. RESULTS AND DISCUSSION 



The calculations of the MD for the various s-p and s-d shell nuclei, with N = Z, have been carried out on the basis 
of Eq. ( |36| ) and the analytical expressions of the one- and two-body terms which were given in Sec. III. Two cases 
have been examined, named case 1 and case 2 corresponding to the analytical calculations with HO orbitals without 
and with SRC, respectively. 

The parameters b and j3 of the model in case 1 and for 4 He, 16 0, 36 Ar and 40 Ca in case 2 were the ones which have 
been determined in our previous work |3l| ] by fit of the theoretical F c h(q), derived with the same cluster expansion, to 
the experimental one. These values of the parameters are given in Table I. The values of the correlation parameter (3 
of the open shell nuclei which have been reported in Ref. |u| were quite large. That is the correlations for these nuclei 
were quite small. The MD of the open shell nuclei, which we found with these values of the parameters, had a high 
momentum tail at values of k larger than expected. As that seems to us quite unreasonable we tried to redetermine 
more carefully the parameters of the model by fit of the theoretical F c h{q) to the experimental one in order to obtain 
a better fit. 

The new values of b and f3 for case 2 and for 12 C, 24 Mg, 28 Si and 32 S are shown in Table I. The theoretical F c h(q) 
for these nuclei, which are shown in Fig. 1, are closer to the experimental data than they were in Ref. 0. From 
the values of x 2 , which have been found in cases 1 and 2 and also from Fig. 1 it can been seen that the inclusion of 
SRC's improves the fit of the form factor of the above mentioned nuclei. Also, all the diffraction minima, even the 
third one which seems to exist in the experimental data of 24 Mg, 28 Si and 32 S are reproduced in the correct place. 

Although the values of the parameters b and (3, for the open shell nuclei, are different from those reported in Ref. 
pl[ , their behaviour, still, indicates that there should be a shell effect in the case of closed shell nuclei. This behaviour 
has an effect on the MD of nuclei as it is seen from Fig. 2, where the MD, of the various s-p and s-d shell nuclei 
calculated with the values of b and (3 of Table I for case 2, have been plotted. It is seen that the inclusion of SRC's 
increases considerably the high momentum component of n(k), for all nuclei we have considered. Also, while the 
general structure of the high momentum component of the MD for A = 4, 12, 16, 24, 28, 32, 36, 40, is almost the 
same, in agreement with other studies [pp| JTl|^0] | , there is an A dependence of n(k) both at small values of k and in 
the region 2 fm _1 < k < 5 fin. . The A dependence of the high momentum component of n(k) is larger in the open 
shell nuclei than in the closed shell nuclei. It is seen that the high momentum component is almost the same for the 
closed shell nuclei 4 He, 16 O and 40 Ca as expected from other studies j|,|||l(|. 

In the previous analysis, the nuclei 24 Mg, 28 Si and 32 S were treated as Id shell nuclei, that is, the occupation 
probability of the 2s state was taken to be zero. The formalism of the present work has the advantage that the 
occupation probabilities of the various states can be treated as free parameters in the fitting procedure of F c h{q). 
Thus, the analysis can be made with more free parameters. For that reason we considered case 2* in which the 
occupation probability 772s of the nuclei 24 Mg, 28 Si and 32 S was taken to be a free parameter together with the 
parameters b and (3. We found that the \ 2 values become better, compared to those of case 2 and the A dependence 
of the parameter (3 is not so large as it was before. The new values of b and (3 are shown in Table I and the theoretical 
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F c h(q) in Fig 1. The values of the occupation probability r]2 S of the above-mentioned three nuclei are 0.19982, 0.17988 
and 0.50921, respectively, while the corresponding values of r]id, which can be found from the values of r]2s through 
the relation 

Vld = [{Z - 8) - 2t7 2s ]/10, 

are 0.36004, 0.56402 and 0.69816, respectively. The MD of these three nuclei together with the closed shell nuclei 
4 He, 16 O and 40 Ca found in case 2 are shown in Fig. 3. It is seen that the A dependence of the high momentum 
component is now not so large as it was in case 2. As F c h{q) calculated in case 2* is closer to the experimental data 
than in case 2, we might say that this result is in the correct direction, that is the high momentum component of 
the MD of nuclei is almost the same. We would like to mention that experimental data for n(k) are not directly 
measured but are obtained by means of y-scaling analysis |27j| and only for 4 He and 12 C in s-p and s-d shell region. 
We expect that the above conclusion could be corroborated if new experimental data are obtained in the future for 
MD for several nuclei and we carry out a simultaneous fit both to MD and to form factors. 

Finally, in table I we give the one and the two-body terms of the mean kinetic energy, (T) , of the various s-p and 
s-d shell nuclei calculated on the basis of Eq. (|39|), as well as the rms charge radii, (t^/i) which are compared 
with the experimental values. It is seen that the introduction of SRC's (in case 2) increases the mean kinetic energy 
relative to case 1 (((T case 2) — (T case i))/(T case 2)) about 50% in 4 He and 23% in 24 Mg. This relative increase follows 
the fluctuation of the parameter (3. Also the values of the kinetic energy in percents, 100(Tsi,;c)/ (Trotai), as well as 
the ratio < Trotai > /(T#o) follow the fluctuation of the parameter (3. In closed shell nuclei there is an increase of 
the above values by the increasing of mass number. 



V. SUMMARY 



In the present work, general expressions for the correlated OBDM and MD have been found using the factor cluster 
expansion of Clark and co-workers. These expressions can be used for analytical calculations, with HO orbitals and 
in principle for numerical calculations with more realistic orbitals. 

The analytical expressions of the OBDM, MD and mean kinetic energy for the s-p and s-d shell nuclei, which have 
been found, are functions of the HO parameter b, the correlation parameter f3 and the occupation probabilities of the 
various states. These expressions are suitable for the systematic study of the above quantities for the N = Z, s-p and 
s-d shell nuclei and also for the study of the dependence of these quantities on the various parameters. 

It is found that, while the general structure of the MD at high momenta is almost the same for all the nuclei we 
have considered, in agreement with other studies, there is an A dependence on n(k) both at small values of k and 
the high momentum component. The A dependence of the high momentum component becomes quite small if the 
occupation probability of the 2s-state for the nuclei 24 Mg, 28 Si and 32 S is treated as a free parameter in the fitting 
procedure of the charge form factor. 
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TABLE I. The values of the parameters b and /3, of the mean kinetic energy per nucleon, (T) and of the rms charge radii, 
( r ch) 1 ^ 2 i f° r various s-p and s-d shell nuclei, determined by fit to the experimental F c h(q). Case 1 refers to the HO wave 
function without SRC and case 2 when SRC are included. Case 2* is the same as case 2 but with the occupation probability 
of the state 2s taken to be as a free parameter. The experimental rms charge radii are from Ref. [41]. 



Case Nucleus b [fm] /3 [fur 2 ] <T) [Mev] (r 2 ch } 1/2 [fm] 

HO SRC Total Theor. Expt. 



1 4 He 1.4320 15.166 15.166 1.7651 1.676(8) 

2 4 He 1.1732 2.3126 22.594 7.310 29.904 1.6234 

1 12 C 1.6251 17.010 17.010 2.4901 2.471(6) 

2 12 C 1.5190 2.7468 19.469 6.111 25.580 2.4261 

1 16 1.7610 15.044 15.044 2.7377 2.730(25) 

2 16 1.6507 2.4747 17.121 6.493 23.614 2.6802 

1 24 Mg 1.8495 16.162 16.162 3.1170 3.075(15) 

2 24 Mg 1.8103 4.2275 16.870 4.239 21.109 3.0948 
2* 24 Mg 1.7473 2.4992 18.109 6.505 24.614 3.0638 

1 28 Si 1.8941 16.099 16.099 3.2570 3.086(18) 

2 2S Si 1.8236 3.0020 17.369 5.564 22.933 3.2159 
2* 28 Si 1.7774 2.4440 18.283 6.922 25.205 3.1835 

1 32 S 2.0016 14.878 14.878 3.4830 3.248(11) 

2 32 S 1.9368 3.0659 15.891 4.976 20.867 3.4425 
2* 32 S 1.8121 2.6398 18.154 6.761 24.915 3.2822 

1 36 Ar 1.8800 17.273 17.273 3.3270 3.327(15) 

2 36 Ar 1.8007 2.2937 18.827 8.590 27.417 3.3343 

1 40 Ca 1.9453 16.437 16.437 3.4668 3.479(3) 

2 40 Ca 1.8660 2.1127 17.863 8.754 26.617 3.5156 
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q [fm" 1 ] q [fm" 1 ] 

FIG. 1. The charge form factor of the nuclei: 12 C (a), 24 Mg (b), 28 Si (c) and 32 S (d) for various cases. Case HO+SRC* 
corresponds to the case when the occupation probability rj2 a is treated as a free parameter. The experimental points of 12 C are 
from Ref. [42] and for the other nuclei from Ref. [43]. 
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1 2 3,4 5 

k [fm ] 

FIG. 2. The correlated momentum distribution for various s-p and s-d shell nuclei calculated with the parameters b and 
of case 2 when the nuclei 24 Mg, 28 Si, 32 S and 36 Ar were treated as lei shell nuclei. The normalization is J n(k)dk = 1. 




1 2 3 4 5 

k [fm 1 ] 

FIG. 3. The correlated momentum distribution for the closed shell nuclei 4 He, 16 0, 40 Ca calculated as in Fig. 2 and for the 
open shell nuclei 24 Mg, 28 Si, 32 S calculated with the parameters b, (3 and rj2 S of case 2* when they were treated as ld-2s shell 
nuclei. The normalization is as in Fig. 2. 
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